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machines (SSMs). Recently, we showed tliat two QSMs ^Al and A^2 witli ni and 77-2 states, 



^^ ' respectively, are equivalent iff they are (ni + 71-2) ^-equivalent (Theoretical Computer Science 



358 (2006) 65-74). However, using this result to check the equivalence likely needs exponential 



pi 1' expected time. In this paper, we consider the time complexity of deciding the equivalence 

pH . between QSMs and related problems. The main results are as follows: (1) We present a 

^ ' polynomial-time algorithm for deciding the equivalence between QSMs, and, if two QSMs 

^'■ 

are not equivalent, this algorithm will produce an input-output pair with length not more 

k><( \ than (ni +n2)^. (2) We improve the bound for the equivalence between QSMs from (ni + 712)^ 

j^ ■ to n^ + 722 — 1, by employing Moore and Crutchfield's method (Theoretical Computer Science 

237 (2000) 275-306). (3) We give that two MO-lQFAs with ni and ^2 states, respectively, are 

equivalent iff they are (ni + 71-2) ^-equivalent, and further obtain a polynomial-time algorithm 

for deciding the equivalence between two MO-lQFAs. (4) We provide a counterexample 

showing that Koshiba's method to solve the problem of deciding the equivalence between 

MM-lQFAs may be not valid, and thus the problem is left open again. 
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1. Introduction 



Over the past two decades, quantum computing has attracted wide attention in the aca- 
demic community J211I29J . To a certain extent, this was motivated by the exponential speed- 
up of Shor's quantum algorithm for factoring integers in polynomial time ^^ and afterwards 
Grover's algorithm of searching in database of size n with only 0{y/n) accesses |19| . 

Quantum computers — the physical devices complying with the rules of quantum mechan- 
ics were first considered by Benioff [Sj, and then suggested by Feynman ^^. By elaborating 
and formalizing Benioff and Feynman's idea, in 1985, Deutsch [T3j re-examined the Church- 
Turing Principle and defined quantum Turing machines (QTMs). Subsequently, Deutsch J14j 
considered quantum network models. In 1993, Yao |SS1 demonstrated the equivalence between 
QTMs and quantum circuits. Quantum computation from the viewpoint of complexity theory 
was first studied systematically by Bernstein and Vazirani |12j . 

Another kind of simpler models of quantum computation is quantum, finite automata 
(QFAs), that can be thought of as theoretical models of quantum computers with finite 
memory. This kind of computing machines was firstly studied independently by Moore and 
Crutchfield 28 , as well as Kondacs and Watrous I^. Then it was deeply dealt with by 
Ambainis and Freivalds |]P, Brodsky and Pippenger TT, and the other authors (e.g., name 
only a few, [2-4,7-10,18], and for the details we may refer to [SJ)- '^^^ study of QFAs is 
mainly divided into two ways: one is one-way quantum finite automata (IQFAs) whose tape 
heads only move one cell to right at each evolution, and the other is two-way quantum finite 
automata (2QFAs), in which the tape heads are allowed to move towards right or left, or to be 
stationary. (Notably, Amano and Iwama 3^ dealt with an intermediate form called l.SQFAs, 
whose tape heads are allowed to move right or to be stationary, and, particularly, they showed 
that the emptiness problem for this restricted model is undecidable.) Furthermore, by means 
of the measurement times in a computation, IQFAs have two fashions: measure-once IQFAs 
(MO-lQFAs) proposed by Moore and Crutchfield [^, and, measure-many IQFAs (MM- 
IQFAs) studied first by Kondacs and Watrous ^^. 



The characteristics of quantum principles can essentially strengthen the power of some 
models of quantum computing, but the unitarity and linearity of quantum physics also lead 
to some weaknesses. We briefly state some essential differences between the QFAs stated 
above and their classical counterparts by two aspects. One is from their power. The class 
of languages recognized by MM- 1 QFAs with bounded error probabilities is strictly bigger 
than that by MO-lQFAs, but both MO-lQFAs and MM-lQFAs recognize proper subclass of 
regular languages with bounded error probabilities [10,11,26,28]. (Also, the class of languages 
recognized by MM- 1 QFAs with bounded error probabilities is not closed under the binary 
Boolean operations JDEl-) Kondacs and Watrous [22j proved that some 2QFA can recognize 
non-regular language L^q = {a"&"|n > 0} with one-sided error probability in linear time 
(Freivalds ^HI proved that two-way probabilistic finite automata can recognize non-regular 
language L^q with arbitrarily small error, but it requires exponential expected time |23) . As 
it is well-known, classical two-way finite automata can accept only regular languages |24j). 

The other difference is from the viewpoint of decidability. By PAix) we denote the 
probability of the automaton A accepting input string x. Then the four cut-point languages, 
recognized by A with cut-point A £ [0,1], are defined by L^^ = {x : Pa{x) ixi A}, for 
ixiG {<,<,>,>}. When A is an MM-IQFA, Blondel et al. jH] proved that the problems 
of determining whether L^^ (ixiG {<,>}) are empty are decidable, but when \x\£ {<,>}, 
such problems are undecidable. In contrast, when A is a probabilistic automaton, all these 
emptiness problems for \x\£ {<,<,>,>} are undecidable [30] . 

Recently, another finding is concerning an essential difference between quantum sequential 
machines (QSMs) and stochastic sequential machines (SSMs). SSMs j^ may be viewed as 
a generalization of probabilistic automata [32,30], since an SSM that has only one output 
element and some accepting states are assigned, reduces to a probabilistic automaton. Two 
SSMs, say A4i and M.2 having ni and n2 states, respectively, and the same input and output 
alphabets, are called equivalent if they have equal accepting probability for any input-output 
pair {u,v). As was known, a crucial result concerning SSMs by Paz [50] is that Mi and M2 
are equivalent iff they are (ni + n2 — l)~equivalent (that is, their accepting probabilities are 
equal for any input-output pair whose length is not more than ni -|-n2 — 1). Recently, Gudder 
j20j first defined sequential quantum machines (SQMs), a quantum analogue of SSMs, and 
Gudder asked whether or not such an equivalence consequence also holds for SQMs. Then, 



Qiu ST re-defined an equivalent version called quantum sequential machines (QSMs), that 
were formally a quantum counterpart of SSMs, just as quantum finite autom,ata (QFAs) to 
probabilistic automata. Qiu [^ demonstrated that there are two QSMs with ni and 71-2 
states, respectively, such that they are (ni + n2 — 1) -equivalent, but not equivalent, after all. 
Hence, Gudder's problem was given a negative answer. 

Latterly, we pTJ further proved that two QSMs Mi and M2 having ni and n2 states, 
respectively, and the same input and output alphabets / and O, are equivalent iff they are 
ijii + 712) ^-equivalent, a new feature in contrast to the (ni + 722 — 1) -equivalence for SSMs 
|3nj . However, using this result to check the equivalence between QSMs needs exponential 
expected time {0{m^'^'^~^^'^' )) where ?n, = |/| x \0\. 

The remainder of the paper is organized as follows. In Section 2, we recall the definitions of 
SSMs and QSMs, and related results. Section 3 is the main part. In Subsection 3.1, we detail 
a polynomial-time algorithm (0(m.(ni + 71-2)^^)) for deciding the equivalence between QSMs. 
In Subsection 3.2, we improve the bound for the equivalence between QSMs from (ni + 77-2)^ 
to n\ + n^ — 1, by employing Moore and Crutchfield's method J2H]- Section 4 is concerning 
the equivalence between one-way QFAs. In Subsection 4.1, we provide a polynomial-time 
algorithm for deciding the equivalence between MO-lQFAs. In Subsection 4.2, we provide 
a counterexample showing that the method stated in ^S] to decide the equivalence between 
MM-lQFAs is not valid. Finally, some remarks are made in Section 5 to conclude this paper. 

2. Preliminaries 

In this section, we briefly review some definitions and related properties that will be used 
in the sequel. 

Firstly, we explain some notations. An n-dimensional row vector (oi 02 . . . a„) is called 
stochastic if Oj > (i = 1, 2, . . . , n), and X^ILi ^« ~ 1' ^^ particular, it is called a degenerate 
stochastic vector if it has 1 only in one entry and else Os. A matrix is called stochastic if its 
each row is a stochastic vector. As usual, for non-empty set /, by |/| we mean the cardinality 
of /, and by /* we mean the set of all finite length strings over /. For u £ I* , \u\ denotes the 
length of u; when \u\ = 0, n is an empty string, denoted by e. We denote I^ = I* — {e}. For 



input alphabet / and output alphabet O, the set of all input-output pairs is defined as 

{{u,v) el* xO* : \u\ = \v\}. 
For any input-output pair (n, f), we denote by l{u,v) the length of u or v. 

For the details on stochastic sequential machines (SSMs), we can refer to [30]. Next, we 
recall the definition of quantum sequential machines (QSMs), a quantum counterpart of SSMs 
[30]. 

Definition 1 ([31]). A QSM is a 5-tuple M = {S , rji^^ , I , O , {A{y\x) : x £ I,y £ O}) where 
S = {si,S2, ■ ■ ■ ,Sn} is a finite set of internal states; r]^^ is an n— dimensional degenerate 
stochastic row vector; / and O are input and output alphabets, respectively; A{y\x) is an 
n X n complex matrix satisfying "^ A{y\x)A{y\x)^ = I for any x £ I, where the symbol f 
denotes Hermitian conjugate operation and I is unit matrix. In particular, for input-output 
pair (e, e), ^(eje) = I. 

Remark 1. It is worth pointing out that, before the definition of QSMs [31], Gudder [20] 
first defined sequential quantum machines (SQMs), an equivalent version of QSMs. The 
equivalence between QSMs and SQMs was proved by Qiu [31]. The reader can refer to 
[20,31] for more information. 

In brief, we may denote QSM Ai as {S, rj^^,!, O, {A{y\x)}). In the QSM Ai defined above, 
if matrix A{y\x) = [aij{y\x)], then aij{y\x) (resp. \aij {y\x)\'^) represents the amplitude (resp. 
the probability) of the machine entering state Sj and yielding y after x being inputted with 
the present state Sj. Thus, given the QSM Ai above, we let Pjii{v\u) denote the probability 
of Ai printing the word v after having been fed with the word u, and it is defined as follows: 

2 

Pm{v\u) = rii^-,A{v\u) , (1) 

where t/jq is the initiation-state distribution of At, (n, v) = (3;iX2 . . . Xm, yiy2 ■ ■ ■ Um) denotes 
an input-output pair, and A{v\u) = A{yi\xi)A{y2\x2) ■ ■ ■ A{ym\xm)- Clearly, we have 

PMiv\u) = r]i,A{v\u)A{v\u)^r]l. (2) 

If the initiation-state distribution 77J,, in QSM At is omitted, then we call A^ uniniti- 
ated QSM (UQSM). For a UQSM A^, by P]^' {v\u) we mean the probability that, with the 
initiation-state distribution rji^ being specified, A4 prints v after u being inputted. 
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Definition 2. Let A^ be a UQSM. Any two initiation-state distributions 77J,, and 77^^, of Ai 
are said to be equivalent (resp. /c-equivalent) with respect to M-, if Pj^'{v\u) = Pj(^^ {v\u) 
for any input-output pair {u,v) (resp. for any input-output pair {u,v) with l{u,v) < k). 

In the following, we define the equivalence between machines. 

Definition 3. Two machines (SSMs, SQMs, or QSMs) Mi and A^2 with the same input 
and output alphabets are called equivalent (resp. ^-equivalent) if Pjii-^{v\u) = PM2i'^W) for 
any input-output pair {u,v) (resp. for any input-output pair {u,v) with l{u,v) < k). 

Remark 2. Given a QSM Ai, from Eq. (1), we know Pjii{e\e) = 1. Therefore, in Definition 
3, we can require l{u,v) > 1. In what follows, for verifying the equivalence between QSMs, 
we only consider the input-output pair {u,v) with l{u,v) > 1. 

A crucial result concerning SSMs is that two SSMs with rii and n2 states, respectively, 
are equivalent iff they are (rii + 722 — l)-equivalent [30]. Therefore, Gudder asked whether it 
holds for SQMs. Then Qiu [31] gave its negative answer. Recently, we further showed that 
two QSMs M-i and M2 are equivalent iff they are (ni -|- 712) ^-equivalent [27]. The two results 
are stated in the following theorem. 

Theorem 1 ([31,27]). (1) There exist QSMs (or SQMs) Mi and M2 with ni and n2 
states, respectively, and the same input and output alphabets, such that though Mi and M2 
are (rii + n2 — 1)- equivalent, they are not equivalent. 

(2) Two machines (SQMs or QSMs) Mi and M2 with ni and n2 states, respectively, and 
the same input and output alphabets, are equivalent iff they are {ui + n2)'^ -equivalent. 

3. Equivalence between QSMs 

In this section, we consider further the equivalence between QSMs. In Subsection 3.1, 
based on the way stated in [27], we give a polynomial-time algorithm that takes as input two 
QSMs and determines whether they are equivalent. In Subsection 3.2, by employing Moore 
and Crutchfield's method [28], we give a better bound for the equivalence between QSMs. 



3.1. A polynomial-time algorithm for the equivalence between QSMs 

As stated before, directly testing Theorem 1 (2) for the equivalence between QSMs needs 
exponential expected time. Therefore, in this subsection, we present a polynomial-time algo- 
rithm for the equivalence between QSMs. 

Before presenting the algorithm, we recall the definition of direct sum of two matrices. 
Suppose that Amn and Bm are m x n and k x / matrices, respectively. Then the direct sum 

A-mn © Bki is an {m + k) x (n + /) matrix, defined as Amn © B^i 



-^mn U 

Bki 



Now, we can present the main theorem as follows. 



Theorem 2. There is a polynomial-time algorithm, that takes as input two QSMs Mi and M2 
and determines whether ^Al and M2 o'^e equivalent. Furthermore, if the two QSMs are not 
equivalent, then the algorithm outputs an input-output pair {u,v) satisfying that Pa4i(w|^) 7^ 
P_\4^{v\u), and l{u,v) < (rii + n2)^, where ui and 722 are the numbers of states in A^i and 
M2, respectively. 

Proof. Given two QSMs having the same input and output alphabets: 

Ml = {Si,r]i^^,I,0,{Ai{y\x)}) and M2 = {S2,Vjo,I,0,{A2{y\x)}), 

where |5i| = ni and |S'2| = n2. We construct UQSM M = {S , I , O , {A{y\x) : x £ I,y £ O}), 
where S = SiU S2, A{y\x) = Ai{y\x) © A2{y\x). For any input-output pair {u,v), denote 

D{v\u) = A{v\u)A{v\u)'< . (3) 

Then for any input-output pairs (n, f ) and {x,y) where l{x,y) = 1, we have 

D{yv\xu) = A{yv\xu)A[yv\xuy 

= A{y\x)A{v\u)(A{y\x)A{v\u))^ 

= A{y\x)D{v\u)A{y\x)^ . (4) 

Let 

D = {D{v\u) : (n, v) is input — output pair, and l{v, u) > 1}, (5) 

and p = (77jy, 0) and p = (0, f/jo), where p and p are (ni + n2)-dimensional row vectors, and 
can serve as two different initiation-state distributions of M . Then for any input-output pair 
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{u,v), we have 



and similarly, 



P^{v\u) = pA{v\u)A{v\u)^ p^ 

= ri^-,Ai{v\u)Ai{v\u)''7]l^ 

= PmiHu), (6) 

P^ivlu) = PmMu). (7) 



Therefore, A^i and M2 are equivalent (i.e., PMiiv\u) = PM2i'v\u) for any input-output pair 
{u,v)) if and only if p and p are equivalent with regard to M, that is, for any D{v\u) £ P, 

pD{v\u)p'< = pD{v\u)p\ (8) 

Let $(2?) be the linear subspace spanned by D, and let ;B be a basis for the subspace ^(V). 
Clearly, the total space as to ^{T>) consists of all (rii + n2)— order complex square matrices, 
together with the usual operations of matrices, whose dimension is (rii + 722)^. Hence B has 
at most (ni + 722)^ elements, and the two QSMs Mi and M2 are equivalent if and only if 
Eq. (8) holds for every vector D{y\u) G B. 

Design of the algorithm. Without loss of generality, we assume that / = {a} and O = 
{0, 1}. Then we define binary tree T as follows. Tree T has a corresponding node D{v\u) 
(defined in Eq. (3)) for every input-output pair {u,v) £ I* x O* . The root of T is D{e\e) 
that is identity matrix /. In the definition of D (Eq. (5)), we notice that l{v,u) > 1, so, we 
exclude -D(e|e) when searching for the basis B of subspace ^(D). However, we still make it 
act as the root of tree T for the sake of convenience. Every node D{v\u) in T has two children 
D{Ov\au) and D{lv\au). For any x £ I, and y £ O, D{yv\xu) can be calculated from its 
parent D{v\u) by Eq. (4). 

Our algorithm is described in Figure 1 which is to efficiently search for the basis B of $(2?) 
by pruning tree T. In the algorithm, queue denotes a queue, and B, the basis stated above, 
is initially set to be the empty set. We visit tree T by breadth-first order. At each node 
D{v\u), we verify whether it is linearly independent of B. If it is, we add it to B. Otherwise, 
we prune the subtree rooted at D{v\u). We stop traversing tree T after every node in T has 
been either visited or pruned. The vectors in the resulting set B will form a basis for ^{"D), 



which will be proven later on. At the end of the algorithm, we verify whether Eq. (8) holds 
for every vector in B. If yes, then the two QSMs are equivalent. Otherwise, the algorithm 
returns an input-output pair {u,v) satisfying that PA4^(f|ii) / Pa/(2(^I^)- 

Figure 1. Algorithm for the equivalence between QSMs. 



Input: Ml = (5i,r7.„,{a},{0,l},{Ai(y|x)}), and M2 = 


= ('S'2,?7jo 


{a},{0,l},{A2{y\i 


c)}) 


Set B to be the empty set; 










queue «— D{€\e); 










while queue is not empty do 










begin take an element D{v\u) from 


queue; 








it D{v\u) ^ span{B) then 










begin add vector D{v\u) to B; 


//D{e\e) is not added to B. 




add D{Qv\au) and D{lv\au 


) to queue; 








end; 










end; 










ifVSe^, (r;,„,0)i?(77,„,0)t = (0,r;, 


g)B{0,r]jgy then return (yes) 




else return (the pair {u,v): {rjig,0)B{r]ig,Oy 7^ (0 


V,o)B{0, 


v^oH 





Remark 3. The basic idea regarding our algorithm is to efficiently search for the basis B of 
$(2?). The foundation of our algorithm is the breadth-first search for traversing a tree. The 
method used in our algorithm is to prune tree T. By pruning some unwanted subtrees, we 
do not need to visit all nodes with height not more than (rii + 712)^ (here the height of root 
is defined as 0), such that we can greatly reduce the number of nodes to be visited. 

Validity of the algorithm. Now we explain why the resulting set B form a basis for ^(D). 
From the analysis above, we know that after running the algorithm, we will get a pruned 
tree. We denote the resulting tree by Tp which is formed by the nodes in the following set 

BU {D{aov\aiu) : D(v\u) £ B, D{aov\aiu) £ span{B),ai £ I,ao£ O}, 

where the former part B consists of the internal nodes of tree Tp, and the latter part comprises 
the leaf nodes of tree Tp. 



For i > 0, we let 

Bi = {D{yv\xu) : D{v\u) is a leaf in Tp, l{x,y) = i}, (9) 

where when i > 1, set Bi consists of unvisited nodes in tree T which have distance i from a 
leaf in tree Tp; when i = 0, set Bq is the set of leaves of Tp. Then it can be readily seen that 

oo 

V = Bu\jBi. (10) 

j=0 
Proving that B forms a basis for ^{T>) amounts to showing that span{B) = span{T>). Equiv- 
alently, we only need to prove the following proposition. 

Proposition 3. For all i> 0, Bi Q span{B). 

Proof. Let B = {i?i, ^2, . . . , 5^} for some m < (rii + 722)^. We show the proposition 
by induction on i. The basic case Bq C span{B) follows straightforward from our analysis 
above. Now assume that Bi C span{B). Then for any input-output pairs {u,v),{x,y), and, 
((To,cJ/), where l{ao,crj) = 1, such that D{v\u) is a leaf and l{x,y) = i, by Eq. (9) we 
know D{yv\xu) £ Bi, and, with the assumption, D{yv\xu) = 'YlT=i'^j^j ^^^ some Qj G C 
(j = 1, 2, . . . , m); furthermore, for any Bj G B, say Bj = D{vj\uj) for some input-output pair 
{uj,Vj), we have A{ao\cri)BjA{ao\<Jiy = D{aoVj\ajUj) G span{BuBo). Therefore, we get 
that 

D{aoyv\aixu) = A{ao\ai)D{yv\xu)A{ao\cri)'' 

= A{ao\fJi)(^ajBj]A{ao\ai)^ 

= J2'^j[M(^oWi)BjA{ao\ai)A 

G span{B U Bq) = span{B). 
This shows that Sj+i C span{B), and the proposition is proved. D 



Complexity of the algorithm. Firstly we assume that all the inputs consist of complex 
numbers whose real and imaginary parts are rational numbers and that each arithmetic 
operation on rational numbers can be done in constant time. Because the basis B has at most 
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(ni + 722)^ elements, the nodes to be visited will be at most 0((ni + 712)"^). (Here we need 
recall a result that to verify whether a set of n-dimensional vectors is linearly independent 
needs time 0{n'^) [17].) At every visited node D{v\u) the algorithm may do two things: (i) 
verifying whether or not the (ni + 71,2) ^-dimensional vector D(v\u) is linearly independent of 
the set B, which needs time 0{{ni + n2)®) according to the result in [17] just stated above; 
(ii) calculating its children nodes by Eq. (4) (if D(v\u) ^ B), which can be done in time 
0((ni + n2)^). Thus the total runtime is 0((ni + n2)^^). 

So far, we have completed the proof of Theorem 2. D 



Remark 4. In the algorithm above, for convince we consider only the case where |/| = 1 
and |0| = 2. In general, let m = |/| x |0|. Then the algorithm almost keeps on except that 
the total nodes to visit will be at most 0{m.{ni +712)^), and as a result, the time complexity 
will be 0{m.{ni + 712)^^). 

3.2. An improved bound for the equivalence between QSMs 

In this subsection, we give an improved bound for the equivalence between QSMs, using 
the bilinearization technique given by Moore and Crutchfield [28]. Firstly, we define a new 
model as follows. 

Definition 4. A bilinear machine (BLM) is a four-tuple M = {S,TT,{M{a)}„^Y.,ri) over 
alphabet S, where S with [S"! = n is a finite state set, tt G C^^", r/ G C"^i and M{a) G C"^" 
for o" G S. 

Associated to a BLM M, the word function fjn : S* ^ C is defined in the form: fjvii'w) = 
■kM{wi) . . . M{wn)ri, where tf = wi . . . ti;„ G S*. 

A probabilistic automaton (PA) is a BLM with the restriction that vr is a stochastic vector, 
r] consists of O's and I's only, and the matrices M{a) (a G S) are stochastic. Then, the word 
function fjn associated to PA M has domain in [0, 1]. 

Definition 5. Two BLMs (include PAs) Mi and M2 over the same alphabet S are said to 
be equivalent (resp. /c-equivalent) if fMii'^) = fM2i'^) ^^ ^^Y w £ T,* (resp. for any input 
string w with |tf | < k). 
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As stated in Paz [30], the result with regard to the equivalence between SSMs can also be 
applied to PAs. Therefore, based on Paz [30], Tzeng [34] considered further the equivalence 
between PAs, giving a polynomial-time algorithm to the problem. Now, the results are stated 
in the following. 

Theorem 4 ([30,34]). Two PAs Aii and M2 with rii and n2 states, respectively, are equiv- 
alent if and only if they are (ni + n2 — 1) -equivalent. Furthermore, there is a polynomial-time 
algorithm that takes as input two PAs Mi and M2 o.nd determines whether Mi and M2 are 
equivalent. 

Remark 5. In fact, one can readily find that Paz 's way [30] can also be applied to BLMs, 
and the algorithm given by Tzeng [34] still works for BLMs. Therefore, Theorem 4 holds for 
the more general model — BLMs. 

Next, we transform a QSM to a BLM by the way given by Moore and Crutchfield [28], 
and then obtain an improved result for the equivalence between QSMs. That is the following 
theorem. 

Theorem 5. Two QSMs Mi and M2 with ni and n2 states, respectively, are equivalent if 
and only if they are (nf + 722 — 1) -equivalent. 

Proof. Given an n-state QSM M = {S, rj^^, I, O, {A{y\x)}), let hj (j = 1, . . . n) be a column 
vector that has only 1 in the jth element and else Os. Then we have 



'2 

Pm{v\u) = r]ioA{v\u) =^\vioMv\u)h 



Oj\ 



i=i 



5^(%. ^ vViMvH ® Aiv\uy] {hj h*) 



i=i 



= {mo ® O [A{v\u) ® A{v\u)*] Y,{hj ® h* 
Here we can construct a BLM M = (S ,iT,M{a)aeT.,'n) ^^ follows: 



S = {{y\x) ■.ye0,xel}, and M{{y\x)) = A{y\x) A{y\x); 
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Then we get PMiui ■ ■■Vmlxi . ■ ■ Xm) = fj^'i{yi\xi) . . . {ym\xm))- Hence every n-state QSM 
can be transformed to an equivalent n^-state BLM, and by Remark 5, we have proved the 
theorem. D 

Remark 6. Considering the equivalence between two QSMs, we have improved the bound 
from (ni + n2)^ to (n^+n^ — 1) by the way of Moore and Crutchfield [28], which seems to imply 
that the way used in Subsection 3.1 is unwanted. Nevertheless, the way used in Subsection 
3.1 offers us a different insight to QSMs and even other quantum computing models, and 
maybe can be used to solve some new problems concerning quantum computing models. 

4. Equivalence between one-way QFAs 

In this section, we consider the equivalence between one-way QFAs. More specifically, 
in Subsection 4.1, we present a polynomial-time algorithm for the equivalence between MO- 
IQFAs, by means of the idea in Subsection 3.1; in Subsection 4.2, we provide a counterexample 
showing that the method used in [25] to decide the equivalence between MM-lQFAs may be 
not valid. 

4.1 Equivalence between MO-lQFAs 

First, we review the definition of MO-lQFAs [28,11]. 

An MO-IQFA ^ is a 5-tuple A = (Q, S, qo, {A{x) : x G S}, F), where Q is a finite set of 
states (let \Q\ = n); qq is the initial state; S is a finite set of input symbols; A{x) denotes 
an n X n unitary evolution matrix for each x £ T,; F <^ Q is the set of accepting states, 
with corresponding projection matrix Pace = diag(po Pi ■ ■ -Pn-i) where for i = 0, . . . , n — 1, 
Pi equals to 1 ii qi £ F else 0. 

As usual, let {qi\ denote the n-dimensional row vector (0 • • • 1 • • • 0) whose {i + l)th com- 
ponent is 1 and the others Os (i = 0, 1, . . . , n — 1). Any configuration of A is described by a 
unit row vector in the superposition form {ip\ = X^^^Tq ai{qi\, with that X^^^Tq |ajp = 1, and 
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Ui denoting the amplitude of A being in state qi. HA is in configuration (^| and reads an 
input symbol a G S, then the new configuration of A becomes {ifj \ = {ip\A{(T). 

The probability of A accepting input string u = X1X2 ■ ■ ■ Xm is defined as 

P^»(^) = IK90|A(u)Paccf (11) 

where A{u) = A{xi)A{x2) ■ ■ ■ A{xm)- For simplicity, we often write P^°(n) by Pj[{u) if no 
confusion results. 

We introduce two definitions regarding the equivalence between MO-lQFAs as follows. 

Definition 6. Two MO-lQFAs Ai and A2 having the same set of input symbols are called 
equivalent (resp. /c-equivalent) if for any input string u (resp. for any input string u with 
\u\ < k), they have equal accepting probability, i.e., PAiiu) = -P42(^)- 

Definition 7. Given an MO-lQFAs A whose initial state is not specified, then two states qi 
and (72 in A are called equivalent (resp. /c-equivalent) if for any input string u (resp. for any 
input string u with \u\ < k), P^(u) = P^{u). 

Concerning the equivalence between MO-lQFAs, we have the following theorem. 

Theorem 6. Two MO-lQFAs Ai and A2 are equivalent if and only if they are (rii + ^2)^- 
equivalent, where ni and n2 are the numbers of states in Ai and A2, respectively. Further- 
more, there is a polynomial-time algorithm that takes as input two MO-lQFAs Ai and A2 
and determines whether Ai and A2 are equivalent. 

Proof. We give a brief proof by four steps below. 

1. Let MO-IQFA A = {Q, S, qo, {A{x) : x G S}, F). Then for any u G S*, we have 

^T(^) = ||(<?0|A(u)Paccf 

= {qo\A{u)P^,,Pl^,A{u)^\qo) 
= {qo\A{u)P^^^A{u)^\qo), 

where notation |-) denotes the conjugate transpose of (-j. 
Denote F{u) = A{u)Pi,ccA{u)'> . Then 

P^°(n) = {qo\F{u)\qo), (12) 

F{xu) = A{x)F{u)A{x)'<. (13) 
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2. Denote 

^ = {F(n) :nG S*}, (14) 

T{k) = {F{u):ue^*,\u\<k}. (15) 

As we did in the proof of Theorem 1, a set of hnearly independent vectors can be found 
in ^(n^) (n = \Q\) such that any vector in ^ is a hnearly combination of these vectors. 
Therefore, by Eq. (12), we obtain that two initial states qq and q^ for A are equivalent 
iff they are n^-equivalent. 

3. As in the proof of Theorem 2, given two MO-lQFAs ^i = (Qi, S, qq, {Ai{x)}, Fi) and 
A2 = {Q2,T.,po,{A2{x)},F2),weletA= {QiUQ2,T.,{Ai{x)®A2{x) : x G S},FiUF2) 
(we assume Qi r\ Q2 = 0). Then the equivalence between Ai and A2 amounts to the 
equivalence between the initial states p and p with regard to A, where p and p have 
corresponding vectors {{qo\,0) and (0, {po\), respectively. 

4. In virtue of the above considerations and the idea in Subsection 3.1, we describe an 
algorithm in Figure 2. Analogous to Theorem 2 and Remark 4, the time-complexity of 
this algorithm is 0{m.{ni + n2)"'^^), where ni = \Qi\, n2 = \Q2\ and m = |S|. 





Figure 2. Algorithm for the equivalence between MO 


-IQFAs. 


Input: Ai = (Qi, go, {0, 1}, {Ai{x) : x = 0, 1} 


Fi), and 




A2 


= (Q2,Po, {0, 1}, {A2ix) : X = 0, 1}, F2) 

Set B to be the empty set; 

queue <— node{e); 

while queue is not empty do 

begin take an element F{u) from queue; 

if F{u) ^ span{B) then 
begin add vector F{u) to B; 

add F{Ou) and F{lu) to queue; 
end; 
end; 








if VB e B, i{qo\,0)B{\qo),0) = (0, (pol)B(O, |po)) then return (yes) 




else return (the string u: {{qo] , 0) B {\qo) , 


0)^(0,(po|)S(0, 


bo))); 
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From the above four steps we have completed the proof. D 

Remark 7. It is worth indicating that Brodsky and Pippenger [11] and Koshiba [25] also con- 
sidered the equivalence problem concerning MO-lQFAs. Their methods can be described by 
two steps: (i) firstly using the bilinearization technique [28] to convert MO-lQFAs to gener- 
alized stochastic finite automata [28] ; (ii) secondly determining the equivalence of generalized 
stochastic finite automata. Their difference is regarding step (ii): Koshiba [25] applied the 
tree pruning technique [34] to determine the generalized stochastic systems' equivalence, while 
Brodsky and Pippenger [11] employed Paz [30] 's method to do that. Therefore, Koshiba [25] 
gave a polynomial-time algorithm for the problem, but Brodsky and Pippenger [11] did not 
consider its efficiency. One can find that our method is different from [25, 11]. 

4.2 A counterexample for the equivalence between MM-lQFAs 

Gruska [22] proposed as an open problem that is it decidable whether two MM-lQFAs are 
equivalent. Then Koshiba [25] tried to solve the problem. For any MM-IQFA, Koshiba [25] 
wanted to construct an equivalent MO-glQFA (like MO-IQFA but with evolution matrices 
not necessarily unitary) and then decided the equivalence between MO-glQFAs using the 
known way on MO-lQFAs. Nevertheless, we find that the construction technique stated in 
[25, Theorem 3] may be not valid, and as a result, the problem is in fact not solved there. We 
will give a counterexample to show its invalidity. In the following, we adopt the definitions 
of QFAs stated in [11] where only the right end-marker symbol $ is considered. So the reader 
can refer to [11] for the definitions and we do not detail them here. 

First let us recall the method stated in [25, Theorem 3] for constructing MO-glQFAs from 
MM-lQFAs. Given an MM-IQFA M = {Q,^,{Ua}aeJ2u{$},Qo,Qacc,Qrej), an MO-glQFA 
Ai = {Q , S, {C/^}o-g2u{$})<?0)-^) is constructed as follows: 

• Q' = Q U {g,, : a G S U {$}}\Q„cc, and F = {g, : a G S U {$}}; 

• U'„\q) = h ai\qi) h aA\qa) when Ua\q) = h ai\qi) h o^I^'a) and qA G Qacc ; 

• add the rules: U^\qa) = {qa) for all [go-) £ -^• 
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Koshiba [25] deemed that the construction technique stated above can ensure that for any 
input word, the accepting probabihty in M is preserved in Ai , which wih be shown to be 
not so. 

Now we turn to the following counterexample provided by us, showing the invalidity of 
the above method. 

A counterexample Let MM-IQFA M = (Q, S, {t/'o-}o-esu{$})'?0) Qacc, Qrej), where Q = 
{QOj Qii Qacc, Qrej} with the set of accepting states Qacc = {Qacc} and the set of rejecting states 
Qrej = {Qrej}', S = {a}; qq is the initial state; {f^cr}cresu{$} are described below. 



Ua{\qo)) = ^\qo) + ^\qi) + ^\ 



iacc/ } 






Uai\qi)) = Trko) 7=\qi) + TilQacc), 



U${\qo)) = \qacc), U${\qi)) = \qrej)- 
Next, we show how this automaton works on the input word aa$. 

1. The automaton starts in |go). Then [/„ is applied, giving ^\qo) H — Tf ki) + ^\qacc)- 
This state is observed with two possible results. With probability (i)^, the accepting 
state is observed and then the computation terminates. Otherwise, a non-halting state 
i|go) + -7^\qi) (unnormalized) is observed and then the computation continues. 

2. After the second a is fed, the state ^\qo) + -75 ki) is mapped to ^(^ + -y^)\qo) + ^(^ ~ 
4=)|gi) + ^(^ + -j^)\qacc)- This is observed with two possible results. With probability 
[2(2 H — 75)]^' t^6 computation terminates in the accepting state qacc- Otherwise, the 

computation continues with a new no-halting state ^(^ + -y^)\qo) + -j^i^ 7n)\qi) 

(unnormalized) . 

3. After the last symbol $ is fed, the automaton's state turns to |(| + -^)\qacc) + ^(| ~ 
~y=)\qrej)- This is observed. The computation terminates in the accepting state \qacc) 
with probability [^(i + 4s)]^ or in the rejecting state \qrej) with probability [-tsC^ — 

The total accepting probability is (i)^ + [l(i + -j^)]"^ + [1(1 + ■j=)]'^ = | + _I^. 
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Now according to the construction technique [25, Theorem 3] stated before, we get an 
MO-glQFA M' = (<9',S,{[/^}^gsu{$}i9o,-^) where Q' = {qo,qi,qrej,qa,q$}, F = {qa,q$} 
and {f^crlo-esulS} ^^e described below. 

U'ailqo)) = -ko) + -/|l9i) + 2^^^)' 

^$(ko)) = \q$), u'^i\qi)) = \qrej), 

U'ailQa)) = \qa), f^$(ka)) = ka)- 

When the input word is aa$, the automaton works as follows. Starting from state \qo), 
when the first a is fed, the automaton turns to state ^\qo) + -y^lqi) + ^\qa) ■ After the second a 
is fed, the state is mapped to ^{^ + ^)\qo) + ^i^-4^)\qi) + [^ + ^{^ + 4^)]\qa)- After the last 
symbol $ is fed, the state is mapped to ^(i + -^)\q$) + ^(^ - ^)krej) + [l + Uh + Tf)]'^"^' 

The total accepting probability is[^(^ + ^)]^ + [^ + |(^ + ^)]^ = | + ^- 

Now it turns out that the accepting probability in the original MM-IQFA is not preserved 
in the constructed machine as expected in [25]. Therefore, the invalidity of the method [25, 
Theorem 3] has been shown. 

Remark 8. (1) The essential reason for the invalidity of the way in [25] is that the accepting 
state set F in Ai does not cumulate the accepting probabilities in the original MM-IQFA. 
Instead, it accumulates just the accepting amplitudes. In addition, we know that in general, 
|ap + |6p 7^ |a + 6p. Therefore, the way in [25] leads to invalidity. (2) Due to the complex 
accepting behavior of MM-lQFAs, it is likely no longer valid to decide the equivalence between 
MM-lQFAs as we did for MO-lQFAs. To our knowledge, so far there seems to be no existing 
valid solution to this problem. Therefore, the equivalence between MM-lQFAs is worth 
considering further. 

5. Concluding remarks 

In this paper, based on the results in [27, 31], we presented a polynomial-time algorithm 
{0{m.{ni -1-^2)"'^^)) for determining the equivalence between two QSMs with rii and n2 states, 
respectively, and, if they are not equivalent, this algorithm will produce an input-output pair 
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with length not more than (ni+n2)^- Furthermore, by using the way of Moore and Crutchfield 
[28], we obtained that two QSMs Ali and A^2 are equivalent iff they are {nf + 71,2 — 1)- 
equivalent, which improves the result in [27]. 

We also proved that two MO-lQFAs ^1 and A2 that have ni and n2 states, respectively, 
and the same input alphabet S with |S| = m, are equivalent if, and only if they are (?^l+?1-2)^~ 
equivalent. In terms of the idea of the algorithm for QSMs, we further provided a polynomial- 
time algorithm {0{m.{ni + n2)^^)) for the equivalence between ^1 and A2- 

In addition, considering the problem of deciding the equivalence between MM-lQFAs, we 
provided a counterexample showing that the method stated in [25] to solve the problem may 
be not valid, and therefore the problem is left open again. 

The further problems are regarding the minimization of states for QSMs [20,31]. As well, 
the equivalence concerning 2QFAs [26] is worthy of consideration. 
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